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The barrier formation for metal/organic semiconductor interfaces is analyzed within the Induced 
Density of Interface States (IDIS) model. Using weak chemisorption theory, we calculate the induced 
density of states in the organic energy gap and show that it is high enough to control the barrier 
formation. We calculate the Charge Neutrality Levels of several organic molecules (PTCDA, PTCBI 
and CBP) and the interface Fermi level for their contact with a Au(lll) surface. We find an excellent 
agreement with the experimental evidence and conclude that the barrier formation is due to the 
charge transfer between the metal and the states induced in the organic energy gap. 



I. INTRODUCTION 



The field of electronic materials based on molecular films is developing very fast. Designing new organic-based 
materials and devices requires a detailed knowledge of basic processes, such as those controlling the formation of 
metal/organic interface barriers 

The evolution of our understanding of this aspect of organic interfaces follows a path which reminds us of the slow 
process of understanding inorganic semiconductor/metal interfaces. Twenty-five years ago, the main problem was to 
understand the mechanism of the formation of the Schottky barrier. Starting with the Schottky and Bardeen models, 
research developed new ideas based on the Defect model and the Induced Density of Interface States model Q . The 
present consensus on inorganic/metal interfaces is that, unless the interface has many defects, the Schottky barrier 
formation is controlled by "intrinsic" interface states induced by the interaction between the inorganic semiconductor 
and the metal. 

Several models have been advanced to explain organic semiconductor/metal interface barrier formation. The 
Schottky-Mott model was originally believed to hold for these interfaces, assuming that no interface dipole is formed 
at the junction, a situation which was subsequently disproved in most cases 0,0- At reactive interfaces, the metal- 
molecule chemical reaction creates gap states that pin the Fermi level, a situation that is analagous to that described 
by the Unified Defect model proposed for inorganic semiconductor/metal interfaces 0. Compression of the metal 
, surface electronic tail by the organic molecules, leading to a change in the metal workfunction, has also been suggested 
f^) ' as a mechanism operating in these interfaces Q- 

This communication focuses on non-reactive interfaces between metals and thin films of low weight organic 
molecules. We study several organic semiconductor/metal interfaces within the IDIS model Q and show that the 
induced densities of states at these junctions are large enough to control the interface barrier formation. We analyze 
the chemical interaction between Au and several organic molecules: 3, 4, 9, 10- perylenetetracarboxylic dianhydride 
(PTCDA), 3, 4, 9, 10- perylenetetracarboxylic bisbenzimidazole (PTCBI) and 4, 4', -N, N'- dicarbazolyl biphenyl 
, (CBP) (Figure^a)). Our quantum-mechanical analysis shows how the weak chemical interaction creates, neverthe- 
less, a local density of states in the organic energy gap, which is large enough to make the IDIS model applicable and 
the definition of a Charge Neutrality Level (CNL) for the organic molecule meaningful. Our results for these junctions 
O ■ explain their pinning behaviour, which is characterized by the interface slope parameter, S = A. low value of S 

, (as is the case of PTCDA/ Au and PTCBI/ Au) corresponds to strong Fermi level pinning, whereas a higher value of 
S (CBP/Au is an example) means that the change of the barrier height with the metal workfunction is larger: as will 
be discussed below, this is associated with a smaller density of states induced in the organic energy gap. 
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II. MODEL AND THEORETICAL SOLUTION 



Consider, initially, the case of a PTCDA/ Au interface. Experimental evidence indicates that PTCDA molecules 
lie flat on the Au surface, and that the PTCDA monolayer has the two-dimensional pattern shown in Figure Hfb). 
PTCDA crystals can be thought of being formed by the repetition of this layer along the direction perpendicular to 
the surface 9j. It is important to realize that, as is typical of organic crystals, intermolecular bonds are weak van 
der Waals bonds that preserve the individuality of the molecules. This simplifies the analysis of the PTCDA-Au 
interaction, reducing it to the case of a single molecule deposited with its plane parallel to the surface. 
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FIG. 1: (a) Organic molecules studied in this paper: PTCDA, PTCBI and CBP; (b) Two-dimensional pattern of a PTCDA 
crystal (from |8|). 

We should stress that, in principle, the molecule-molecule interaction induces some (small) broadening of the 
electronic levels of each individual molecule, but does not create an electron density of states in the molecular energy 
gap. Since the Schottky barrier formation depends on these gap states, we can safely neglect the molecule-molecule 
interaction and consider only the single molecule-metal interaction. 

Our analysis of the organic semiconductor/metal interface is made in several steps. First, the organic molecule is 
analyzed using a DFT-LCAO method [Hj ; this is a DFT-based theory, which uses a local-orbital basis and the orbital 
occupation for describing the exchange and correlation energies of the system. It has been shown elsewhere that this 
approach is equivalent to other more conventional DFT methods that use a local exchange and correlation energy 

El- 

The main problem with DFT methods for organic molecules (and other small molecules) is that their one-electron 
eigenvalues do not represent real electron or hole excitations. In particular, the DFT energy gap is not directly related 
to either the transport or the optical gap. The advantage of our approach is that, by using a variation of Koopman's 
theorem [Tl| , we can calculate the molecular electronic levels by introducing a relaxation energy in the DFT method 
that is directly related to the exchange-correlation energy. We have shown for PTCDA IFJ t hat, although this effect 
is important, the relaxation of the molecular orbitals calculated within DFT is negligible jlJj. This means that, after 
the DFT-LCAO calculation for the organic molecule, we need to renormalize the energy levels by means of Koopman's 
relaxation energies in order to obtain a realistic energy spectrum. We can keep, however, the electronic wavefunctions 
of the molecule obtained from the DFT calculation. 

Figure |21 illustrates this situation for PTCDA. Panel (a) shows the energy spectrum calculated within DFT-LCAO, 
while panel (b) shows the spectrum after introducing Koopman's relaxation energies. Note how the PTCDA energy 
gap increases from 1.6 eV to 5.7 eV due to these relaxation effects; the molecular wavefunctions associated with the 
ionization (HOMO) and affinity (LUMO) levels are assumed to be practically the same before and after the energy 
relaxation. 

There are also solid-state effects associated with long-range electronic polarization. These are mainly associated 
with screening of the metal and the organic crystal of the electronic field created by the extra charge (electron or 
hole) introduced in the molecule |l3j |. This correction reduces the energy gap by ~1.5eV, although the ionization 
and affinity wavefunctions of the molecule are not expected to present important modifications. Other polarization 
effects, due to lattice relaxation and vibronic effects only introduce small corrections, around 0.2eV, further reducing 
the molecular energy gap [||. Regarding PTCDA, we fit the molecular transport gap to 3.2eV, which presumably 
takes into account all the effects discussed above Q; this is done by introducing a rigid shift between the empty and 
occupied states of Figure Hfb). For PTCBI and CBP, we fit the transport gaps to 3.1 eV and 5.1 eV respectively, 
values which were deduced from experimental optical gaps and exciton binding energies [l4j |. 

In a second step, we calculate the induced density of states at the organic semiconductor/metal interface using 
chemisorption theory in the limit of weak interaction between the two systems [l6j . In our model, we assume the 
organic molecule (PTCDA, PTCBI or CBP) to be deposited flat on the metal surface, at a distance d from the last 
metal layer, which we will take for the rest of this paper to be Au(lll). 

In our analysis, we start with the organic molecule wavefunctions, tpi, obtained from the DFT-LCAO method 
discussed above, and the metal density of states matrix, p a p(E), where a and (3 refer to the local-orbital basis used 
to describe the metal properties, which are calculated using the DFT local-orbital code Fireball ^5|. In the limit of 
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FIG. 2: DFT-LCAO energy spectrum (long bars: it states; short bars: a states) before (a) and after (b) Koopman's relaxation 
energies; the relaxation in the wavefunctions, however, is small. The DOS is obtained by introducing a 0.5 eV FWHM Gaussian 
broadening for each level. 

weak PTCDA-metal interaction, the main effect of the metal is to broaden the molecular levels Ei by the quantity I\ 

T l = 2TrJ2 \T lv \ 2 SiE, - Ei), (1) 

where Ti V is the hopping interaction between the molecular orbital ipi and the metal eigenfunction, ?/;„. Equation 
^ can be rewritten in a more convenient way by using the local-orbital basis for the molecule and the metal. Writing 
tpi = Cij<f>j, equation Stakes the form 

r. ( = 2ir c ij Tj a PapiEi) Tpji Cj'i. (2) 

jj'a/3 

Neglecting off-diagonal terms with j ^ j' and a ^ ft, which tend to cancel each other out, equation is further 
simplified to 

F, = 27T \Ci,j\ 2 \Tj,a\ 2 pa, a ( E i)' ( 3 ) 

In our calculation, we have only included the interaction of the Au 6s orbital with the different orbitals of the 
organic molecules, C 2s and 2p, N 2s and 2p, O 2s and 2p, and H Is. This implies that in equation^ a refers only to 
the Au 6s orbitals. Figure [3] shows these interactions as a function of the Au-atom distance, which obviously depends 
on the organic-metal distance, d. 

The organic-metal separation is a difficult issue. First, no experimental value exists for these systems. Second, 
conventional DFT codes cannot be used to calculate d, due to the weak van der Waals interaction between the metal 
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FIG. 3: Hopping elements between the atoms which make up the organic molecule and Au as a function of distance. 



CNL (eV) El(eV) 



PTCDA 


2.5 


3.2 


PTCBI 


2.4 


3.1 


CBP 


3.0 


5.1 



TABLE I: Charge Neutrality Levels (measured from the center of the HOMO) and peak-to-peak transport gaps, for the three 
organic materials. 



and the organic molecule 9]. Using indirect information based on: (a) the PTCDA-PTCDA stacking distance, d ~ 3.2 
A, and, (b) the atomic radius of Au (~ 0.5 A larger than that of C), we assume the distances between the last Au 
layer and the plane of the organic molecules considered in this paper to be around 3.5 ± 0.3 A. We make use of this 
value to calculate Ti from equation 

Once Ti is calculated, each molecular level Ei contributes to the organic LDOS with the Lorentzian function 

I T J1 (4) 

The total LDOS is obtained by adding all the contributions from the different molecular orbitals. 



III. RESULTS AND CONCLUSIONS 



Figure H shows our results for PTCDA, PTCBI and CBP. In these calculations d = 3.5 A, rj ~ 0.5 eV and 
Tf ~ 0.25 eV. Note the different broadening of the it and a levels. These values were obtained using equation |3 and 
the hopping integrals of Figure El for d = 3.5 A. Note also the uncertainty on the broadening, due to the uncertainty 
in the metal-organic distance (see below). 

Figure 2] shows the different molecular levels for each molecule, the 7r-((7-)states being drawn as long (short) bars. 
For each molecule, the CNL is calculated by imposing charge neutrality conditions: the total electronic charge below 
the CNL integrates to the number of occupied molecular states. Table U shows the calculated CNLs (measured from 
the ionization level). Also included in this table are the transport gaps for the sake of comparison. Note the similarity 
between PTCDA and PTCBI, with the CNL level rather close to the affinity level. This is due to the similarity of 
their energy spectra and transport gaps. In both cases the CNL is close to the affinity level because of the distribution 
of 7r-states around the energy gap. The large density of states below the ionization level pushes the CNL upper in 
the gap. For CBP, the CNL is closer to the molecular midgap, due to the larger energy gap and the more symmetric 
distribution of 7r-states around the HOMO and LUMO. 
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FIG. 4: (TOP to BOTTOM:) IDIS, CNL and interface Fermi level for PTCDA/Au, PTCBI/Au and CBP/Au. Long (short) 
bars correspond to the it (a) states neglecting the metal-molecule interaction. The insets show the region around Ef in detail. 



The interface slope parameter, S — jx^i is given by: 

S dEF 1 (5) 

d<p M l + ATre 2 D{E F )5/A 1 y ' 

where D(Ep) is the induced density of states at the Fermi energy (twice the values shown in Figure 0}, d is the 
metal-organic distance and A is the area associated with one organic molecule (see Figure ^ A=120 A 2 , 191 A 2 and 
251 A 2 , for PTCDA, PTCBI and CBP respectively). Table ILTI compares the experimental and calculated values for 
the slope parameter. Note again the similarity between PTCDA and PTCBI, and the larger value found for CBP. 
This is due to the smaller D(Ep) and the larger area A. The different pinning behaviour of the three interfaces is 
explained: the low value of S, as is the case of PTCDA and PTCBI on Au(lll), corresponds to a high pinning at the 
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PTCDA/Au PTCBI/Au CBP/Au 

S (theory) 0.2 0.2 0.5 
S (exp) 0.0 00 0.6 

TABLE II: Calculated and experimental values for the interface slope parameter, S. 





PTCDA/Au 


PTCBI/Au 


CBP/Au 


E F (theory) (eV) 


2.5 


2.3 


2.5 


E F (exp) 


2.5 


2.1 


2.4 



TABLE III: Theoretical and experimental interface Fermi level positions, measured from the center of the HOMO. 

organic CNL. For CBP/Au(lll), on the other hand, the larger value of S reflects the higher ability of Ep to move 
within the organic energy gap. 

This can be seen in the relation CNL — Ep = S{GNL — 4>m)- Having calculated CNL and S, we obtain the interface 
Fermi level straightforwardly. The pinning at the interface reduces the initial difference CNL — (j)M, to the injection 
barrier Ep — (f>M- 

Figure0]also shows the position of the Fermi level for interfaces with Au. The position, measured from the ionization 
level, is shown for each molecular film in Table ITTT1 The agreement between the theoretical and experimental positions 
of Ep and values of S is remarkable, although a small difference appears for the slope parameter. This presumably 
reflects the approximations introduced in our calculation. The main source of inaccuracy comes from the assumed 
value of d, which has an error bar of around ±0.3 A. This inaccuracy is mainly reflected in the calculated values of 
D(Ep), while the CNL is probably very insensitive to that modification. To ascertain this point, we recalculated the 
CNL changing V* and Tf by factors of up to 2, and found that its value remains practically the same. This is not 
the case for D(Ep), which changes by 50% when d changes ±0.2 A. Our results suggest that d has probably been 
overestimated for PTCDA and PTCBI by at least 0.3 A in our calculations. 

We stress, however, that our results for the interface Fermi level show very good agreement with the experimental 
data. Moreover, the main outcome of our analysis is to show that the induced density of interface states is large 
enough to play a crucial role in the formation of the metal/organic semiconductor barriers. This allows us to conclude 
that the mechanism associated with the formation of these interface barriers is the charge transfer between the two 
materials due to the weak chemical interaction: this creates an electrostatic interface dipole which tends to align the 
metal Fermi level and the organic CNL. 
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